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We present a new generalized topological current in terms of the order parameter field to 
describe the topological defect system in 0{n) symmetric time-dependent Ginzburg-Landau model. 
With the aid of the ^-mapping method, the structure of the topological defects and the topological 
quantization of their topological charges in TDGL model are obtained under the condition that the 
Jacobian J(^) 7^ 0. We show that the topological defects are generated from the zero points of 
the order parameter field (j), and the topological charges of these topological defects are topological 
quantized in terms of the Hopf indices and Brouwer degrees of 0- mapping under the condition. When 
J(^) — 0, it is shown that there exist the crucial case of branch process. Based on the implicit 
function theorem and the Taylor expansion, we detail the bifurcation of generalized topological 
' current and find different directions of the bifurcation. The topological defects in TDGL model are 

found splitting or merging at the degenerate point of field function <j) but the total charge of the 
topological defects is still unchanged. 
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CN . I. INTRODUCTION 

> 

Q,^ ■ The importance of the role of time-dependent Ginzburg-Landau (TDGL) model in understanding a variety of 
' problems in physics is clear The study of the growth kinetics of systems subjected to rapid temperature quenches 
CO recently been extended to include systems with more complex order parameter symmetries [HJ^ . In particular 

there has been progress on the study of the n-dimensional vector model with nonconserved Langevin dynamics. An 
interface description and numerical simulations of a TDGL equation were used to investigate the intra-surface kinetics 
of phase ordering on toroidal and corrugated surfaces 0. Two-dimensional XY models with resistively-shunted 
I junction (RSJ) dynamics and TDGL dynamics were simulated and it was verified that the vortex responses can be 
^ ■ well described by the Minnhagen phenomenology for both types of dynamics Q . 

G I In recent years, much work has been done on the topological defects in the TDGL model p|,P-pl]|. In certain 
: cosmological |l2| and phase ordering |]l3| problems key questions involve an understanding of the evolution and 
(-H ■ correlation among defects like vortices, monopoles, domain walls, etc. In studying such objects in field theory questions 
O , arise as to how one can define quantities like the density of topological defects and an associated velocity field. Liu 
O ■ and Mazenko have discussed the problem, but unfortunately, for the lack of a powerful method, the topological 
^ ] properties of these systems are not very clear, some important topological information has been lost, also the unified 
theory of describing the topological properties of all defect system in TDGL model is not established yet. 

In our previous work Jl^ , we discussed how one could use the (^-mapping topological current theory to study the 
^ . topological structure of the point defects in phase-ordering systems. We deduced a topological current of point-defect 
■ " ' system in terms of the order parameter field in the context of a d-dimensional 0{n) symmetric TDGL model, where 
d is the spatial dimensionality. Using the expression of the topological current, for point defects (n = d), the point- 
defect velocity field was identified, the topological structure and the topological charges of the point-defect system 
were studied, also the branch process of the point defects was discussed systematically. This analysis will be extended 
here to the case of arbitrary dimensional defect system in TDGL model where n ^ d — k, i.e. to study the case for all 
n < d. 

In this paper, in the light of 0-mapping topological current theory p5[ |, a useful method which plays a important 
role in studying the topological invariants 1 16 1^ and the topological structures of physical systems 1 18 1^ , we will 



investigate the topological quantization and the branch process of arbitrary dimensional topological defects in 0{n) 
symmetric TDGL model. We will show that the topological defects are generated from the zero points of the order 
parameter field 0, and their topological charges are quantized in terms of the Hopf indices and Brouwer degrees of 
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(^-mapping under the condition that the zero points of field (j> are regular points. While at the critical points of the 
order parameter field 0, i.e. the limit points and bifurcation points, there exist branch processes, the topological 
current of defect bifurcates and the topological defects split or merge at such point, this means that the topological 
defects system is unstable at these points. 

This paper is organized as follows. In section 2, we investigate the topological structure of the defect system in 
TDGL model and give the expression of the defect velocity field. In section 3, we will point out that the topological 
charges of these defects are the Winding numbers which are determined by the Hopf indices and the Brouwer degrees 
of the (/)-mapping. In section 4, we study the branch process of the defect topological current at the limit points, 
bifurcation points and higher degenerated points systematically by virtue of the </)-mapping theory and the implicit 
function theorem. 



II. TOPOLOGICAL STRUCTURE OF THE DEFECT SYSTEM IN TDGL MODEL 

We consider a TDGL model for an rt-component order parameter = (0^ (r, t), (/'"(r, t)) in d-spatial dimensions 
{d ^ n — k) governed by the Langevin equation 



dt 



5F 

= K = -T—+ri (1) 



where F is a kinetic coefficient and ^ is a thermal noise which is related to F by fluctuation-dissipation theorem. F 
is a Ginzburg-Landau effective free energy assumed to be of the form 

F^J d''r[^c{W$'f + V{$)] (2) 

where c > and the potential V is assumed to be of the degenerate double- well form. This model is to be supplemented 
by random, uncorrelated, initial conditions. We assume that there is a rapid temperature quench from a high 
temperature to zero temperature where the noise fj in can be set to zero. In the scalar case {n = 1) such 
system order through the growth of domains separated by sharp walls. As time evolves these domains coarsen and 
order grows to progressively longer length scales. In the case of systems with continuous symmetry (n > 1) the 



disordering elements |20 2^ will depend on n and spatial dimensionality d. Thus, for example, for n = d one has 
point defects (vortices and monopoles),for n = d — 1 one has vortex line or string-like defects while for n = d — fc one 
has fc-dimensional topological defect objects. For n > d there are no stable singular topological objects. 

As is well known, the n-component order parameter fleld t) determines the defect properties of the system, and 
it can be looked upon as a smooth mapping between the {d + l)-dimensional space-time X and the n-dimensional 
Euclidean space i?" a.s (p : X ^ i?". By analogy with the discussion in our previous work , from this 

(/)-mapping, one can deduce a topological tensor current as 



1 

A(5'"-i)(^ - 1) 



/^o, = 0, 1, ai, ...,a„ = 1, (3) 

to describe the fc-dimensional topological defects in TDGL model, and the density tensor of the topological defect 
system is defined as p'^i "^* = jOMi -/'fe_ jj^ ^jjjg expression, stands for d/dx^^, ^(5""^) = 27r"/^/F(n/2) is the area 

of (n — l)-dimensional unit sphere S'^~^ and n°'{x) is the direction field of the n-component order parameter field (p 

= \\m\\ = V¥{xmx) (4) 

with n"-{x)n°-{x) — 1. It is obviously that n°-{x) is a section of the sphere bundle S{X) and it can be looked upon 
as a map of X onto a {n — l)-dimensional unit sphere 5'"^^ in order parameter space. Clearly, the zero points of 
the order parameter field 0(x) are just the singular points of the unit vector n°'{x). It is easy to see that j^" "^': are 
completely antisymmetric, and from the formulas above, we conclude that there exists a conservative equation of the 
topological tensor current in (||) 

a^^jMo-M. ^ 0, i = 0,...,fc. 
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and from which we have 



which is just the continuity equation satisfied by p^^'"^'' . 

In the foUowing, we will investigate the intrinsic structure of the generalized topological current j^oA'i-- A»fc i^y making 
use of the ^-mapping method. From , we have 

1 19 1^" 



which should be looked upon as generalized functions . Due to these expressions the generalized topological current 
can be rewritten as 



^ J A(S"-i)(n-2)(ri-l)!' 



2ir ' 



n > 2 
n = 2 



where Cm is a constant 



and G„(||0||) is a Green function 



Defining general Jacobians J^^oAii - Mfc following 

and by making use of the n-dimensional Laplacian Green function relation in 0-space ||l5| ] 

where = { g^ag^a ) is the n-dimensional Laplacian operator in 0-space, we do obtain the (5-function structure of 
the defect topological current rigorously 

This expression involves the total defect information of the system in TDGL model and it indicates that all the defects 
are located at the zero points of the order parameter field (j){x). It must be pointed out that, comparing to similar 
expressions in other papers, the results in is gotten theoretically in a natural way. From this expression, the 
density tensor of jt^ofj-i-'-t^k jg a,lso changed into a compact form 



X 



where D^'i - a^'c (|) = jOMi - Pk(|), We find that jMoMi -p*. ^ q, p^i '^fe ^ only when = 0, which is just the 
singularity of jf^ofj.i---f^k ^ detail, the Kernel of the (/)-mapping is the singularities of the topological tensor current 
j/io/ii - Mfc in X, i.e. the inner structure of the topological tensor current is labeled by the zeroes of ^-mapping. We 
think that this is the essential of the topological tensor current theory and (/)-mapping is the key to study this theory. 

From the above discussions, we see that the kernel of (/)-mapping plays an important role in the topological tensor 
current theory, so we are focused on the zero points of (j) and will search for the solutions of the equations <f>°'{x) = 
(a = 1, ...,n) by means of the implicit function theorem. These points are topological singularities in the orientation 
of the order parameter field 4'{x). 
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Suppose that the vector field (j){x) possesses I zeroes, according to the implicit function theorem [Q, when the 
zeroes are regular points of 0-mapping at which the rank of the Jacobian matrix [dp^cf)""] is n, the solutions of = 
can be expressed parametcrizedly by 

xf" ^ z'^{t,u\- ■ ■,u''), (7) 

where the subscript i represents the i-th solution and the parameters (/ = l,...,fc), combining with the time 
parameter t, span a (fc + l)-dimensional subspace which is called the i-th singular subspace Ni in the space-time X 
corresponding to the (/)-mapping. These singular subspaces Ni are just the world volumes of the topological defects, 
the parameters play the role of the spatial parameters of the topological defects. For each singular subspace Ni, 
we can define a normal subspace Mi in X which is spanned by the parameters {A — 1,- ■ -^n), and the intersection 
point of Mi and Ni is denoted by Pi which can be expressed parametcrizedly by = pf. In fact, in the words of 
differential topology, Mi is transversal to Ni at the point pi. By virtue of the implicit function theorem at the regular 

point Pi, it should be held true that the Jacobian matrices J{^) satisfies J(^) = ^ ^■ 

On the other hand, putting the solutions (0) into (j){x), we have ${zi) = 0, from this expression, and since we 
expect the instantaneous velocity to be orthogonal to the local orientation of the topological defects, we can define 
the velocity via 

that is, the topological current and the density tensor should satisfy 

This expression for the velocity can be used to find the defect velocity distribution in the case of phase-ordering 
kinetics for a non-conscrvcd order parameter p|,p|p4|. It is very useful because it avoid the problem of having to 
specify the positions of the topological defects explicitly. The positions are implicitly determined by the zeros of the 
order parameter field. The general expression with J^^i " '*'' (0/a;) should be useful in looking at the motion of the 
topological defects in TDGL model in the presence of external fields beyond a growth kinetics context 

If we restrict the discussion to the simplest case of point defects {n = d), the corresponding topological current is 
reduced to = S {(f>) J ^ {(f) / x), and Eq.(^ can be put into a conventional form j'' = v^p with v^^ taking the expression 
of v^^ — J^{(l)/x)/ D{4)/x) |14|. Also, the string in TDGl model, which was discussed by Mazenko |^^, is a special 
case for n = d — 1 in our theory. 



III. TOPOLOGICAL QUANTIZATION OF THE DEFECT CHARGES IN TDGL MODEL 

In the following, we will investigate the topological charges of the topological defects and their quantization. Let 
Hi be a neighborhood of pi on Mi with boundary dY^i satisfying pi (jz. dY^i, n S ^ =0. Then the generalized winding 
number Wi of n°'{x) at pi can be defined by the Gauss map n : dT^ S""^^ 

W. = /^^^ n*(e.......n-rf.- A ... A d.-) (10) 

where n* denotes the pull back of map n. The generalized winding numbers is a topological invariant and is also 
called the degree of Gauss map It means that, when the point covers 9Si once, the unit vector will cover 
a region ri[9Si] whose area is Wi times of A(5'"'~^), i.e. the unit vector n°- will cover the unit sphere 5'""^ for Wi 
times. Using the Stokes' theorem in exterior differential form and duplicating the above process, we get the compact 
form of Wi 

W,^ I Si$)Ji^)rv. (11) 

By analogy with the procedure of deducing S{f{x)), since 

r/ T-N _ / +00, for (t){x) = _ / +00, for x £ N , . 

~ I 0, for ${x) ^0^1 0, forxiNr 

we can expand the (5-function 5{(j)) as 
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5(,?) = ^c,<5(iV,), (13) 



where the coefheients must be positive, i.e. Ci =| q |. S{Ni) is the (5-function in space-time X on a submanifold Ni 



d{N,)^ 5''{x~z,{t,v},---,u''))dtd''u. (14) 
Substituting ([l3| ) into ([ll|), and calculating the integral, we get the expression of c,; 

Pi PiTji 



(15) 



where Pi = \Wi\ is a positive integer called the Hopf index pTf of 0-mapping on M^, it means that when the point v 
covers the neighborhood of the zero point pi once, the function (j) covers the corresponding region in 0-space Pi times, 
and = signJ{^)p^ = ±1 is the Brouwer degree of 0-mapping Substituting this expression of q and ( |l3|) into 
(a), we gain the total expansion of the topological current 



PiVi 



or in terms of parameters — (t, w^, • • •, u", , ■ ■ ■, u^) 

From the above equation, we conclude that the inner structure of j^o*'i "^fc or j^o-^i '^fc is labeled by the total 
expansion of 5{<j)), and it just represents I (fc)-dimensional topological defects with topological charges gi = PiTji 
moving in the (d + l)-dimensional space-time X . The (k + 1) -dimensional singular subspaces Ni (i = 1, • • •/) are their 
world sheets in the space-time. Mazenko and Halperin also got similar results for the case of point-like defects 
and line defects, but unfortunately, they did not consider the case Pi 1. In fact, what they lost sight of is just the 
most important topological information for the charge of topological defects. In detail, the Hopf indices Pi characterize 
the absolute values of the topological charges of these defects and the Brouwer degrees rji = +1 correspond to defects 
while = — 1 to antidefects. Furthermore, they did not discuss what will happen when rji is indefinite, which we will 
study in detail in section 4. 

IV. THE BRANCH PROCESS OF THE TOPOLOGICAL CURRENT 

With the discussion mentioned above, we know that the results in the above section are obtained straightly from 
the topological view point under the condition J(0/w)|p. ^ 0, i.e. at the regular points of the order parameter field (f>. 
When the condition fails, i.e. the Brouwer degree rji are indefinite, there should exist some kind of branch processes 
in the topological current of the topological defect system in TDGL model. In what follows, we will study the case 
when J((/>/w)|pi = 0. It often happens when the zero points of field include some branch points, which lead to the 
bifurcation of the topological current. 

In this section, we will discuss the branch processes of these topological defects. In order to simplify our study, let 
the spatial parameters be fixed, i.e. to choose a fixed point on the topological defect. In this case, the Jacobian 

matrices J^°^^"'^''i^) are reduced to 



V y y y V 



A, B = 0, 1, ...,n, = n+ l,...,d, (17) 
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for = {A<n), ip = t, y"^^ = (/ > 1). The branch points are determined by the n + 1 equations 

0"(t,w\---,z;",ti) = O, a=l,...,n (18) 

and 

r+\t,v\---,v\u)^J{^)^Q (19) 

for the fixed u. and they are denoted as {t* ,pi). In 0-niapping theory usually there are two kinds of branch points, 
namely the limit points and bifurcation points ||29|] , satisfying 

^'(|)l(t*.P.)^0 (20) 

and 

■^'(|)I(*-,pO=0: (21) 
respectively. In the following, we assume that the branch points {t*,pi) of 0-mapping have been found. 

A. The branch process at the limit point 

We first discuss the branch process at the limit point satisfying the condition (pO|). In order to use the theorem 
of implicit function to study the branch process of topological defects at the limit point, we use the Jacobian J^{y) 

instead of J{^) to discuss the problem. In fact, this means that we have replaced the parameter thy v^. Then, taking 
account of the condition(|2^) and using the implicit function theorem, we have an unique solution of the equations 
( plsj ) in the neighborhood of the limit point [t* ,Pi) 

t^t{v\u), v' = v\v^,u), i = 2,3, ...,n (22) 

with t* ~ tjpj, u). In order to show the behavior of the defects at the limit points, we will investigate the Taylor 
expansion of (22) in the neighborhood of {t*,pi). In the present case, from ( po|) and (|l9|), we get 



i.e. 

Therefore, the Taylor expansion of (p2|) at the point {t*,pi) gives 



which is a parabola in the — t plane. From (23), we can obtain the two solutions v^-^^^{t,u) and v^2)i^t^)j which 
give the branch solutions of the system ( p^ at the limit point. If j£yr^\{f ,pi) > 0, we have the branch solutions for 

t > t* (Fig 1(a)), otherwise, we have the branch solutions for t < t* (Fig 1(b)). The former is related to the creation 
of defect and antidefect in pair at the limit points, and the latter to the annihilation of the topological defects, since 
the topological current of topological defects is identically conserved, the topological quantum numbers of these two 
generated topological defects must be opposite at the limit point, i.e. + (32V2 = 0. 
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B. The branch process at the bifurcation point 



In the following, let us consider the case (|2l|), in which the restrictions of the system ( |18| ) at the bifurcation point 
(t*,pi) are 

J(-)l(r,p.)=0, J\t)\ ^^^^=0. (24) 
V y 

These two restrictive conditions will lead to an important fact that the dependency relationship between t and is 
not unique in the neighborhood of the bifurcation point {t* ,pi). In fact, we have 



which under the restraint ( |24| ) directly shows that the tangential direction of the integral curve of equation (^5|) is 
indefinite at the point {t* ,pi). Hence, (^ ) does not satisfy the conditions of the existence and uniqueness theorem of 
the solution of a differential equation. This is why the very point (i*, Zi) is called the bifurcation point of the system 
®- 

As we have mentioned above, at the bifurcation point {t*,pi), the rank of the Jacobian matrix [^] is smaller than 
n. For the aim of searching for the different directions of all branch curves at the bifurcation point, we firstly consider 
the rank of the Jacobian matrix [S^l is rt — 1. The case of a more smaller rank will be discussed in next subsection. 
Let Ji(^) = [ipW (a = 1, • • •, n — 1; A = 2, • • •, n) be one of the (n — 1) x (n — 1) submatrix of the Jacobian matrix 
[^] with det Ji(^) 7^ at the point {t*,pi) (otherwise, we have to rearrange the equations of (18)), where (f)'^ stands 
for {dcjf" / dv^). By means of the implicit function theorem we obtain one and only one functional relationship in the 
neighborhood of the bifurcation point [t* .pi) 

^ f^{v\t,u), ^ = 2,3,...,n. (26) 

We denote the partial derivatives as f^^^^ fA^9^^ fA ^ jA^^^^A^d^^ ^^^^^ ^ (||)^ 

we have for a = 1 , . . . , n — 1 

ck'^ = <ly''{v\f{v\t,u),...,r{v\t,u),t,u) = 

which gives 

^d^^' « = 1' 1 (27) 

A=2 



By differentiating (|27| ) and (|2^ ) with respect to and t, and applying the Gaussian elimination method, we can 
find the second order derivatives /j^, /{^ and f^. The above discussions do not matter to the last component 
(j)^{v^, • • •, f ", i, m). In order to find the different values of dv^ /dt at the bifurcation point, let us investigate the Taylor 
expansion of 4>"{v^, ■ ■ u) in the neighborhood of {t*,pi). Substituting ( p^ ) into 0"(u^, • • •,ti",t,u), we get the 

function of two variables and t 

F{t, v\u) = r\v\ f{v\t, u), r{v\t, u), t, u) (29) 

which according to (|l^) must vanish at the bifurcation point 

F{t*,p,)^Q. (30) 



From (29), we can calculate the first order partial derivatives of F{t,v^,u) with respect to and t respectively at 
the bifurcation point it* ,pi) 
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1^ = 0? + E '/'a/A § = c^7 + j:rAff- (31) 



A=2 A=2 



By making use of (|27| ) and (^8|), with the Cramer's rule, the first equation of (24) is expressed as 

l^det Ji(^)|(t.,p,) = 0. 

Since det Ji(^)|(t._p;) 7^ 0, the above equation leads to 

dF 

Q^\it',P.)-0- (32) 
dF 

^l(*%P.)=0. (33) 

The second order partial derivatives of the function F{t,v^,u) are easily to find out from ( |3l| ) which at {t*,pi) are 
denoted by 

d'^F d'^F d'^F 



With the same reasons, we can prove that 



Then, by virtue of (30), (p^), (33) and (p^), the Taylor expansion of F(t,v^,u) in the neighborhood of the bifurcation 
point {t*,pi) gives 

A{v^ - pi)2 + 2B{v^ - pi)(t - r) + C{t -t*f^ 0. (35) 

Dividing ( |35| ) by (w^ — p^)^ or (i — i*)^, and taking the limit t ^ i* as well as p\ respectively, we get two 

equations 



and 

C{^,? + 2B^^+A = Q. (37) 

So we get the different directions of the branch curves at the bifurcation point from the solutions of (^6|) or (|3^). 
There are four possible cases: 

Firstly, A ^ 0, A = 4(5^ - AC) > 0, from Eq. (|3^) we get two different solutions: dv^/dt |i,2= {-B ± 
\/B^~~AC)/A, which is shown in Fig. 2, where two topological defects meet and then depart at the bifurcation 
point. Secondly, A ^ 0, A = A{B^ — AC) = 0, there is only one solution: dv^/dt = —B/A, which includes three 
important cases: (a) two topological defects tangentially collide at the bifurcation point (Fig 3(a)); (b) two topological 
defects merge into one topological defect at the bifurcation point (Fig 3(b)); (c) one topological defect splits into two 
topological defects at the bifurcation point (Fig 3(c)). Thirdly, A = 0, C 7^ 0, A = 4{B^ - AC) > 0, from Eq. (|^) we 
have dt/dv^ — and —2B/C. There are two important cases: (i) One topological defect splits into three topological 
defects at the bifurcation point (Fig 4(a)); (ii) Three topological defects merge into one at the bifurcation point (Fig 
4(b)). Finally, A — C — 0^ Eqs. (^6|) and ( |37| ) give respectively dv^ /dt — and dt/dv^ = 0. This case is obvious as 
in Fig. 5, which is similar to the third situation. 

In order to determine the branches directions of the remainder variables, we will use the relations simply 

dv"^ = ffdv^ + ffdt, A = 2, 3, n 
where the partial derivative coefficients and //^ have given in (^^ and (p8|). Then, respectively 

dv^ „A pA 



dv^ dv^ 

or 



where partial derivative coefficients fi and /j are given by (27) and (|28|). From this relations we find that the values 
of dv^/dt at the bifurcation point {t*,Zi) are also possibly different because (Is^) may give different values of dv^/dt. 
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C. The branch process at the higher degenerated point 



In the following, let us discuss the branch process at a higher degenerated point. In the above subsection, we have 
analyzed the case that the rank of the Jacobian matrix [d(j)/dv\ of the equation ( [l9| ) is n — 1. In this section, we 
consider the case that the rank of the Jacobian matrix is n — 2 (for the case that the rank of the matrix \d(j)/dv\ is 
lower than n — 2, the discussion is in the same way). Let the (n — 2) x (n — 2) submatrix J2{^) of the Jacobian matrix 

[d(j)/dv\ be 



V 



in — 2 in — 2 
»3 



J 



and suppose that det J2{^)\{f .pi) 0- With the same reasons of obtaining (pq), we can have the function relations 

= f'^{v\v'^,t,u), A = 3,4, (39) 



For the partial derivatives ff, /2 and , we can easily derive the system similar to the equations (27) and (Eq), in 



which the three terms at the right hand of can be figured out at the same time. In order to determine the 2-order 
partial derivatives f^, f^, f^, and f^, we can use the method similar to the above mentioned. Substituting 
the relations (|3^) into the last two equations of the system ( p^ , we have the following two equations with respect to 
the arguments w^, f^, t,u 



(40) 



Fi(ui,w2,i,M) = </>"-i(wi,u2,/3(ui,w2,t,u),- • •,/"(w\i;^i,u),i,u) = 
F2iv\v^,t,u) = q^^{v\v^,f{v\v^,t,u),---,riv\v^,t,u),t,u)^0. 

Calculating the partial derivatives of the function Fi and F2 with respect to v^, and i, taking notice of ( ^9|) and 
using six similar expressions to (p2) and (p3), i.e. 



OF, 



0, 



dFj 



dF, 



1,2, 



we have the following forms of Taylor expressions of Fi and F2 in the neighborhood of {t*,pi) 
F,{v\v\t,u) « A,^{v^ - p\f + A,2{v^ - p]){v^ - pi) + A,^{v^ - p\) 



(41) 



*\2 



{t-t*)+ A,^{v' ~pir+ A,5{v' -pf)(t-t*) + A,e{t - n 







J = 1,2. (42) 
In the case of Aji ^ 0, Aj^ ^ 0, by dividing (|2|) by (t — t*Y and taking the limit t t*, we obtain two quadratic 

It and ^ 



equations of ^jr and ^jr- 



~dt 



dv^ dv^ LIU 



dt 



,dv^ 
'dt 



dv^ . 



(43) 



J = 1,2. 

Eliminating the variable dv^ /dt, we obtain a equation of dv"^ /dt in the form of a determinant 



All A12Q + A23 AuQ^ + A15Q + A16 

All A12Q + A13 AiiQ^ + A15Q + Aie 

A21 A22Q + A23 A24Q^ + A25Q + A2e 

A■2^ A22Q + A23 A2iQ^ + A25Q + A26 



A21 J-^22W -^ ^23 

where Q = dv^ /dt, which is a ^th order equation of dv^ /dt 







(44) 
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(45) 



Therefore we get different directions at the higher degenerated point corresponding to different branch curves. The 
number of different branch curves is four at most. If the degree of degeneracy of the matrix [^] is more higher, i.e. 
the rank of the matrix [^] is more lower than the present (n — 2) case, the procedure of deduction will be more 
complicate. In general supposing the rank of the matrix [^] be (n — s), the number of the possible different directions 
of the branch curves is 2^* at most. 

At the end of this section, we conclude that there exist crucial cases of branch processes in our theory of topological 
defect system in TDGL model. This means that a topological defect, at the bifurcation point, may split into several 
(for instance m) topological defects along different branch curves with different charges. Since the topological current 
is a conserved current, the total quantum number of the splitting topological defects must precisely equal to the 
topological charge of the original defect i.e. 



Vi 



for fixed i. This can be looked upon as the topological reason of the defect splitting. Here we should point out that 
such splitting is a stochastic process, the sole restriction of this process is just the conservation of the topological 
charge of the topological defects during this splitting process. Of course, the topological charge of each splitting 
defects is an integer. 
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FIGURES' CAPTIONS 



Fig. 1. (a) The creation of two topological defects, (b) Two topological defects annihilate in collision at the limit 
point. 

Fig. 2. Two topological defects collide with different directions of motion at the bifurcation point. 

Fig. 3. Topological defects have the same direction of motion, (a) Two topological defects tangentially collide at 
the bifurcation point, (b) Two topological defects merge into one topological defect at the bifurcation point, (c) One 
topological defect splits into two topological defects at the bifurcation point. 

Fig. 4. (a) One topological defect splits into three topological defects at the bifurcation point, (b) Three topological 
defects merge into one topological defect at the bifurcation point. 

Fig. 5. This case is similar to Fig. 4. (a) Three topological defects merge into one topological defect at the 
bifurcation point, (b) One topological defect splits into three topological defects at the bifurcation point. 
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